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ABSTRACT 

A dramatic increase in the accuracy and statistics of space-borne cosmic ray (CR) measurements has yielded 
several breakthroughs over the last several years. The most puzzling is the rise in the positron fraction above 
~10GeV over the predictions of the propagation models assuming pure secondary production. The accuracy 
of the antiproton production cross section is critical for astrophysical applications and searches for new physics 
since antiprotons in CRs seem to hold the keys to many puzzles including the origin of those excess positrons. 
However, model calculations of antiproton production in CR interactions with interstellar gas are often employ¬ 
ing parameterizations that are out of date or are using outdated physical concepts. That may lead to an incorrect 
interpretation of antiproton data which could have broad consequences for other areas of astrophysics. In this 
work, we calculate antiproton production in pp-, pA-, and .4.4-interactions using EPOS-LHC and QGSJET-II- 
04, two of the most advanced Monte Carlo (MC) generators tuned to numerous accelerator data including those 
from the Large Hadron Collider (LHC). We show that the antiproton yields obtained with these MC generators 
differ by up to an order of magnitude from yields of parameterizations commonly used in astrophysics. 

Subject headings: cosmic rays 


1. INTRODUCTION 

Antiprotons in cosmic rays (CR) are produced in CR in¬ 
teractions with interstellar gas and are, therefore, called sec¬ 
ondary. The same interactions produce charged and neutral 
mesons that decay to secondary electrons and positrons and 
7 -rays. However, in contrast to CR positrons that can be pro¬ 
duced copiously in pulsars, there is no known astrophysical 
source of primary antiprotons. Since the antiproton “back¬ 
ground” is fairly low this provides an opportunity to search 
for a new phenomenon or an exotic signal if the CR antiproton 
flux is measured accurately. An even smaller background is 
expected for antideuterons, which pr ovide, therefore, another 
target in searches for new physics (Chard onnet et all 1 1 997b 
Dona toet al.ll 200 (I) . 

The first detec tions of antiprotons in C Rs were reported 
35 yea rs ago by Bogom olov et alJ (119791) and iGolden et alJ 
(119791) . Both experiments reported a_ flux that_ w as a 
factor of ~3 above expectations (iGaisser & Levy! 119741) 
and caused a stir at that time even though the reported 
excess was only at the ~3cr level. Now it is clear 
that these meas u remen ts and a followup measurement by 
iBuffington et alj (119811) were likely plagued by a back¬ 
ground of negatively charged particles. Meanwhile, these 
measurements stimulated an intensive discussion about the 
possiblity to observe primary antiprotons produced in ex¬ 
otic processes, such as the annihilation of primordial black 
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CAPRICE98 (Boezi oet al.l2001 ), and most notabl y the BESS 


and BESS-Polar (1993-2008) (lOrito et~akl l2000fc lAbe et alJ 


120121 ) experiments made a set of accurate measurements that 
motivated also new attempts to make more accurate theoreti¬ 


cal evaluations of the flux of secondary antiprotons produced 
in CR interactions with interstellar gas. As a result, it has 
become clear that the measured antiproton flux and p/p ratio 


are generally consistent with the secondary antiproton pro- 

duction (Gaisser & Schaefer 19921 Moskalenko et al. 

1998; 

Simon et al. 1998; Bottino et al. 1998; 

Bieber et al. 

1999; 

Bergstrom et al. 199Donato et al] 2001 

; Moskalenko et al. 

2002 ), but do not rule out the possibility of a weak exotic sig- 

nal. 


More realistic calculations involving the fully numerical 
CR propagation model GALPROP (IStrong & Moskalenkol 
Il998c IStrong et ali 120071) have shown that the pre- 
dicted flux depends on the assumed propagation model 
([Moskalenko et all |2002|) . even though the propagation pa¬ 
rameters were derived from fits to the B/C and 10 Be/ 9 Be 
ratios and the CR fluxes were tuned to local measurements 
(IStrong et al.ll2007l) . In particular, a standard stochastic reac¬ 
celeration model was found to underproduce antiprotons by 
~40% at 3 GeV, while a semi-empirical plain diffusion model 
overproduces them by ~20%. A diffusion-convection model 
was found to be consistent with BESS 95-97 data. The ori¬ 
gin of the remarkable discrepancy with the predictions of the 
physically motivated reacceleration model remained unclear, 
but it was natural to blame the model deficiency rather than to 
declare an exotic signal. One of the possibilities to reconcile 
the model predictions with the data is to assume that some 
fraction of carbon and other primary nuclei at low energies is 
local (iMoskalenko et al.112003b . In this case the appropriate 
decrease in the derived diffusion coefficient increases the pro¬ 
duction of antiprotons, enough to reproduce the data. Another 
solution is to account for the back reaction of CR protons onto 
the interstellar turbulence dPtuskin et alJl2006l) . The stochas¬ 
tic acceleration of CRs by MHD waves is accompanied by 
the so-called damping of the waves, since the wave energy 
is dissipated. This also leads to the effective decrease in the 
diffusion coefficient at a few GV that increases the antiproton 
production. 
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Recent discoveries in astrophysics of CRs, however, have 
changed the landscape dramatically. The most important was 


a clear confirmation of the risi 


PAMELA (lAdriani et ali 12009 


in the positron fraction by 


2 013all that was_earlier no¬ 


ticed by the HEAT collaboration dBeattv et all [2004). This 
was subsequent ly confirmed and extend ed to higher energies 
by Fermi -LAT (lAckermann et all 1201 21) . and measured with 
even greater precision and extended to even higher energies 
by AMS-02 (iAguilar et all 120131) . Most recently, the mea¬ 
surement of the positron fraction has been extended up to 
500 GeV by the AMS-02 collaboration hinting for some flat¬ 
tening of the fraction above 200 GeV (lAccardo et al.ll20l4h . 
but the statistics is still low in this energy range. New 
accurate antiproton measurements were done by PAMELA 
(lAdriani et al]l2009bl 120101 1201 3bl) covering the range from 
70MeV - 200 GeV, and more is expected from AMS-02. 
Above a few GeV, the data are consistent with secondary pro¬ 
duction, in strong contrast with positron results. Another im¬ 
portant milestone is an accurate measurement of the B/C ra¬ 
tio up to lOOGeV/nucleon by PAMELA (lAdriani et alJl2014h 
that is consistent with preliminary AMS-CL0 measurements 
reaching ~400 Ge V/nucleon. Both measurements indicate 
the index of the diffusion coefficient a « 0.4 or even smaller 
(lAdriani et al.ll2014h that supports a Kolmogorov-type power 
spectrum of interstellar turbulence (iKolmogorovll 1941» . thus 
favoring the stochastic reacceleration model for interstellar 
propagation. Besides providing more accurate data over a 
wider energy range these new measurements give an impor¬ 
tant inside into CR acceleration and propagation processes. 
Also relevant, but controversial are the new measurements of 
CR proton and He spectra, the most abundant CR species. 
Combi ned measurement s by PAMELA (lAdriani et al.ll20111) . 
ATIC (iPanov et al.ll2009l) . and CREAM dAhn et al. 120101) hint 
on a break in both spectra at about the same rigidity ~230 
GV. The flattening of the CR proton spectrum at high energies 
was also confirmed by the Fermi -LA T through observations o f 
7 -ray emission of the Earth’s limb (Acker mann et al .1120141) . 
However, preliminary results by AMS-O^do not show any 
feature in the CR proton and He spectra up to ~2 TV. 

Antiproton data and their correct interpretation hold the 
key to the resolution of many astrophysical puzzles. If the 
rise in the positron fraction is due to WIMP annihilations, 
antiprot on data provide important constrai nts on WIMP 
mo dels dD onato et aPl2009l: iCirelli et afll2009l) . for a review, 
see IPorter et all (120m If the rise is due to conventional 
astrophysics, antiproton and B/C measurements extended 
to higher energies may be able to discriminate between 


to ntgner ener gies may be able to d iscriminate between 
the p u lsar (lAronsI 1 198 lE lHarding & RamatvHl987l: iBoularesI 
[19891: iHooner et al.l 120091 iDi Mauro et al l 12071 and SNR 
hypotheses ([Berezhko et all 120031: iBlasil 120091) . The lat¬ 


ter proposes a secondary component with a hard energy 
spectrum that is produced in a SNR shock by accelerated 
protons. It also predicts a rise in all secondaries, _such as 
the v !v and B/C ratios , at hig h energies (Berezhk o et al.l 
Mil IB Iasi & Serpicol 120091 ICholis & Hooperi 120141 

Mertsch & Sarkar 2014). _ However, other a uthors 

(iKachelrieB et al.l 1201 It [KachelrieB & Ostapchenkol 120131) 
pointed out that the results of their time-dependent Monte 
Carlo simulations predict a flattening rather than a distinct 

1 Relative to exp ectati ons as if all positrons in CRs are secondary 
iProtheroe 1982; Moskalenko & Strong 1998). 

2 http://www.ams02.org/2013/07/new-results-from-ams-presented-at-icrc- 
2013/ 


rise. 

The accuracy of the antiproton production cross section is 
critical for astrophysical applications and searches for new 
physics. This is especially true in view of many expec¬ 
tations connected with the upcoming data releases by the 
AMS-02 experiment operating at the International Space Sta¬ 
tion, and by soon-to-be-launched ISS-CREA.V0, CALElQ 
and GA PSQcxpcri merits. This holds even more for new op¬ 
portunities that would open up with antideuteron detection 
in CRs. In turn, the calculation of antideuteron produc¬ 
tion relies on the inclusive antiproton production cross sec- 
tions and the detailed knowledge of two-parti cle correlations 
(iKadastik et alJlMH: IPaT & KachelriefiEoIl) . 

In this work, we analyze antiproton production in pp-, pA-, 
and /l/l-interactions using EPOS-LHC and QGSJET-II-04, 
two of the most advanced Monte Carlo (MC) generators tuned 
to numerous accelerator data including those from the Large 
Hadron Collider (LHC). The antiproton yields obtained with 
these MC generators differ by a factor of few from yields 
of parameterizations commonly used in astrophysics. The 
article is structured as follows: In Section [2] we compare 
the pros and cons of parameterizations and MC generators 
developed specifically for CR interactions at low- and high- 
energies. In Section^ we present a tune of the fragmentation 
procedure used in the QGSJET-II-04 model which leads to an 
improved description of particle production at low energies 
in the presence of high-energy thresholds, as in the case of 
antiprotons. Then we compare in Section [4] the results of the 
modified QGSJET-II-04 and the EPOS-LHC models to exper¬ 
imental da t a and to th e predictions of t he par ameterizations of 
iTan & Ngl (119831) and lDunerrav et all (120031) . Finally, we dis¬ 
cuss the nuclear enhancement of the p yield by nuclear species 
in Section [5]before we conclude. 

2 . P PRODUCTION: MODELS AND 
PARAMETERIZATIONS 

Calculations of secondary CR fluxes, both for astrophysi¬ 
cal applications and searches for new physics, are typically 
based on empirical parameterizations of accelerator data for 
the production spectra of secondary particles in proton-proton 
interactions. Despite the evident convenience of using such 
parameterizations, this practice has a number of caveats. First 
of all, one has to rely on empirical scaling laws when these 
parameterizations are extrapolated outside of the kinematic 
range of the data they are based on. In particular, the high en¬ 
ergy extrapolations of existing parameterizations prove to be 
unreliable, as we will demonstrate in the following. This fail¬ 
ure partly reflects the scarcity of relevant experimental data 
and thus the poor theoretical understanding of the high energy 
behavior of hadronic collisions at the time when the scaling 
laws used as input were developed. Secondly, it is problem¬ 
atic to account for the contributions of nuclear species in CRs 
and/or in the interstellar medium (ISM) using such parame¬ 
terizations: While particle production in proton-nucleus colli¬ 
sions has been studied by a number of fixed target accelerator 
experiments, yielding the spectra of secondaries in the for¬ 
ward hemisphere in the center-of-mass system (c.m.s.), there 
is practically no experimental information on the forward par¬ 
ticle spectra in nucleus-proton and nucleus-nucleus collisions. 
One usually attempts to solve the problem by introducing em- 

3 http://cosmicray.umd.edu/iss-cream/ 

4 http://calet.phys.lsu.edu 

5 http://gammaO.astro.ucla.edu/gaps 
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pirical “nuclear enhancement” factors, which vary, however, 
substantially between different publiciUions. Moreover, as it 
has been demonstrated in IKachelriefi et al.l (120141) . properly 
calculated nuclear enhancement factors depend strongly both 
on the production spectra of the respective particles and on 
the spectral slopes of CR species. 

These difficulties fully apply to the case of antiproton pro¬ 
duction by CRs, which motivated us to reconsider the prob¬ 
lem employing MC generators of hadronic interactions. Com¬ 
paring the results of a number of MC models used in low 
energy nuclear physics to available accelerator data, we ob¬ 
served strong (up to an order of magnitude) deviations from 
the measured p spectra. Hence, we decided to turn to MC 
generators of high energy i nteractions, notably, to those used 


in the high energy CR field (Kalmvkov et al. 

1992 Ahn et al. 

2009; Ostanchcnko 2006; Werner et al. 2006 

). Our choice is 


motivated by the fact that such models have been calibrated 
on numerous accelerator data over a wide energy range and 
survived many consistency tests in high energy CR studies. 
Moreover, most of those models are in a g ood agre ement 
with various data from LHC Run I (d’Enterri a et all201 lh and 
some of them have been r ecently updated ( Ostanchenk ol201 ll 
mil Pierog et alll2013l) . including a re-tuning of model pa¬ 
rameters, based on LHC data. 

However, the application of such MC generators to low en¬ 
ergy hadronic interactions which dominate the production of 
secondary p in the physically important range of kinetic en¬ 
ergies Ep m < lOGeV is rather unwarranted. Such low en¬ 
ergies imply an extrapolation of the underlying theoretical 
approach outside of its range of applicability: Several reac¬ 
tion mechanisms like contributions of Reggeon exchanges, 
intranuclear cascading, etc., which are important at low ener¬ 
gies are irrelevant for high energy interactions and are there¬ 
fore typically neglected. Comparing the model predictions 
with experimentally measured antiproton spectra for proton- 
proton and proton-nucleus collisions for incident momenta 
Piab 100GeV/c, we observed generally a strong disagree¬ 
ment w ith datJT. In the particula r case of the QGSJET-II-04 
model (lQstapchenkoll201H12013 ). this disagreement was es¬ 
pecially large close to the kinematical threshold for p produc¬ 
tion. Only the EPOS-LHC model dPierog et al .11201 3l) demon¬ 
strated a generally reasonable behavior in the low energy 
limit, as will be demonstrated in the following. 

3. LOW ENERGY EXTENSION OF THE QGSJET-II 
MODEL 

The results of EPOS-LHC and QGSJET-II-04 for p produc¬ 
tion are rather similar in the high energy range where both 
MC generators are tuned to LHC data, but deviate strongly in 
the low energy limit. At first sight, this seems surprising as 
their treatment of hadronic interactions at relatively low colli¬ 
sion energies is very similar. The common underlying physics 
includes multiple scattering processes, which are described 
within the Reggeon Field Theory framework (Gribov! 1 196 8) 
by multiple exchanges of Pomerons, i.e. composite states with 
vacuum quantum numbers. Particle production is treated as 
the formation and break-up of strings of color field, which is 
performed using string fragmentation. However, the proce¬ 
dures used in the two models to fragment color strings differ, 
and these differences are responsible for the discussed varia- 

6 Similar discrepancies have been reported i n [Ibarra & Wild <20131 ) for the 
case of the DPMJET-III model fRoesler et al.1200 ll) which is employed in the 
popular FLUKA code l iBattistoni et alJ2007t> . 


tions at low energies. Since QGSJET-II-04 reproduces well 
pion and photon production data down to rather low energies, 
Piab~10GeV/c (IKachelriefi & Qstapchenkoll2012h . it is natu¬ 
ral to suspect that the discrepancies in p production are caused 
by threshold effects related to the relatively high antiproton 
mass. 

These observations motivated us to improve the low energy 
behavior of QGSJET-II, modifying its string fragmentation 
procedure. It is worth stressing that we did not aim at the 
development of a full-scale model for low energy hadronic 
collisions. Instead, we intended a more reliable extrapola¬ 
tion of the current model towards low energies—in order to 
reach an acceptable agreement with experimental data in the 
energy range relevant for calculations of the CR antiproton 
flux. In particular, the modifications introduce no additional 
adjustable parameters. They only modify slightly the string 
fragmentation algorithm, in such a way that the changes have 
no significant influence on the results in the high energy range. 

To describe these modifications in some more detail, 
let us briefly discuss the hadronization procedure used in 
QGSJET-II. The standard physics picture for the string 
break-up is the neutralization of the color field via the 
creation of quark-antiquark and diquark-antidiquark pairs 
from the vacuum, followed by their conversion into 
final-state hadrons (iAndersson et al.l 119831) . In QGSJET- 
II, this_ proc ess is modelled using an iterative proce¬ 
dure (iKalmvkov & Qstanchenkol 19931) . whose parameters are 
determined by the intercepts of the corresponding Regge 
trajectories dKaidalov! 11987 ). After creating a new quark- 
antiquark (diquark-antidiquark) pair, a hadron is formed and 
the reminder of the string proceeds to the next step: Depend¬ 
ing on the string mass, either the same procedure is repeated 
or a two-particle decay is modeled. In the new treatment, 
we introduced an additional weight into the sampling pro¬ 
cess of kinematic variables for each subsequent string break¬ 
up, which is proportional to the two-particle decay phase 
volume evaluated for the mass squared of the reminder of 
the string. While being of minor importance for high mass 
strings, stretched over long rapidity intervals, which are typ¬ 
ically produced in high energy collisions, this modification 
enhances the production of light hadrons (mostly pions) in 
the fragmentation of strings of small masses, at the expense 
of the heavier hadrons. Additionally, the parameters of the 
hadronization procedure, notably the string mass threshold 
for proceeding to a two-particle decay and the relative prob¬ 
abilities for creating quark-antiquark and diquark-antidiquark 
pairs from the vacuum, have been readjusted in order to stay 
in agreement with high energy data. 

In Figs. |T| and [2j we compare the results of the modi¬ 
fied QGSJET-II model (hereafter referred to as QGSJET-IIm) 
and the EPOS-LHC model with selected benchmark accel¬ 
erator datt0. The data shown in Fig. D] are the momen¬ 
tum spectra of p in the laboratory frame in ////-collisions and 
p Be-collisions at _p i a b = 19.2GeV/c (lAllabv et alJ 1197(1 : 
Amaldi et al] 119751) for different angles 0 in the laboratory 
frame. The data shown in Fig. [2] are the Feynman ./.'-spectra 
of antiprotons, l/7ra"E dcr/dxp, where .te = 2 E*/^fs, xf = 
2p*/^/s, with E* and p* being the c.m.s. energy and the ,2- 
component of the momentum, in yy- and y C-collisions a t 
Piab = 158GeV/c (lAnticic et al J l2?)TT)L iBaatar et alJ 120131) . 
In addition, to demonstrate that the introduced modifications 

7 A comparison with additional experimental data sets is presented in Ap¬ 
pendix [a] 
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have not spoiled the treatment of high energy interactions, we 
compare in Fig. [3] the calculated transverse momentum spec¬ 
trum of antiprotons in pip-collisions at \fs = 900 GeV with 
the data of the ALICE experiment (Aamodt et alJl2011 ). As 
can be judged from the Figures, there is a reasonable overall 
description of p-production over a wide energy range. The 
differences with the results of the EPOS-LHC model can be 
used as a measure for model uncertainties in the calculations 
of secondary antiproton spectrum using QCD MC genera¬ 
tors. For comparison, we plot in Figs. [2] and 0 also the re¬ 
sults obtained using the original QGSJET-II-04 model. At 
= 900 GeV, the differences between the models QGSJET- 
Ilm and QGSJET-II-04 are pretty small. On the other hand, at 
lower energies there is a significant reduction of the antipro¬ 
ton yield predicted by QGSJET-IIm, which reaches ~ 20% 
already at 158GeV/c. 


4. Z-FACTORS FOR P PRODUCTION: COMPARISON 
OF MODEL PREDICTIONS 


To compare the impact of different interaction models and 
parameterizations on the predicted CR antiproton spectrum, 
it is c onvenient, similarly to the 7 -ray case (Kachelri efi et alJ 
120141) . to use the corresponding “Z-factors.” They are de¬ 
fined as the spectrally averaged energy fraction transferred to 
antiprotons in proton-proton, proton-nucleus, nucleus-proton, 
or nucleus-nucleus collisions, assuming that the spectra of CR 
species in the relevant energy range can be approximated by a 
power-law, h(E) = Ki E~ ai . Then the contribution qd (Ep) 
to the flux of secondary CRs (here, antiprotons) from inter¬ 
actions of the CR species i with ISM component j ( i,j = 
protons or nuclei) of number density rij, 


Qp ( E p) = n i 


dE 


da i ^P(E,E p ) 

dE e 


Ii(E), ( 1 ) 


can be rewritten as (Kachelriefi et al.f2014l) 


Qp ( e p) = h(E p ) Zp 3 (E p , a,). (2) 


Here, we expressed the Z-factor Zp via the inclusive spectra 
of antiprotons da ( E, zp)/dz p , Zp = Ep/E, as 


Zp (E p , a) = 


a/if). <3) 

t dz 


Note that E corresponds to the energy per nucleon for nuclear 
projectiles. 

The Z-factors Zp depend clearly both on the p production 
spectra and on the spectral slopes on, containing all the depen¬ 
dences on hadronic interaction models. On the other hand, 
these factors are independent of the CR abundances. Using 
two different interaction models. Ml and M2, the ratio of the 
respective contributions to the secondary fluxes equals the ra¬ 
tio of the corresponding Z-factors [c.f. Eq. ©]: 


^/(Ml )(Ep) _ Zp(Ml )(Ep,&i) 

1p(M2) (Ep) ^p(M2) {Ep, Oi) 

In the following, we are going to compare the factors Z p ‘ 
obtained with the modified QGSJET-IIm model and with 
EPOS-LHC to the Z-factors calculated using some widely 
used parameterizations of p-spectra for pp-collisions. How¬ 
ever, before doing so, let us investigate which projectile ener¬ 
gies contribute mainly to Zp P (Ep, a). To this end, we plot 


in Fig. Q] the spectrally-weighted (for definiteness, we use 
a = 2.6) distribution of the energy fraction Zp = Ep/E trans¬ 
ferred to antiprotons 

1 dZ P p P ol _ 1 da i ^(Ep/zp,Zp) 

Zf dz p Zp dz p ( ’ 

for different kinetic energies Ep ln = Ep — m p . For com¬ 
parison, the same distribution for 7 -ray production, 1 / Zl‘/‘ x 
dZP p /dz~f (z 1 = Ery/E), is also shown. Clearly, the range 
of CR proton energies E = Ep/Zp which contributes signifi¬ 
cantly to antiproton production is substantially narrower com¬ 
pared to the 7 -ray case. Since the p spectrum is much softer 
than that of gammas, the region of moderately large z con¬ 
tributes much less to the production of p’s than to the produc¬ 
tion of 7 ’s. Additionally, the distribution becomes substan¬ 
tially narrower with decreasing Ep, which is a consequence 
of both threshold effects and the interaction kinematics: For 
small Ep ln , the region of not too small z p becomes kine¬ 
matically forbidden, because the corresponding proton energy 
E = Ep/zp falls below the production threshold. On the 
other hand, the contribution of the region zp —> 0 is strongly 
suppressed by the spectral factor z p ~ l , c.f. Eq. ©. More¬ 
over, for a given CR proton energy E, antiprotons are pro¬ 
duced most copiously in the c.m.s. central region (xp ~ 0 ), 
which corresponds to Ep ~ in the lab. frame (where 

m^v is the nucleon mass). Thus, the region Ep -C \J2m^E 
or Zp <C \j2rriM/E corresponds to the target fragmentation 
region in the c.m.s., while Ep ln —► 0 corresponds to the kine¬ 
matic boundary (xp —> —1 in c.m.s.), where the p spectrum 
falls down rapidly. This can be clearly seen in Fig. [2] (right) 
for the case of proton-carbon collisions. 

The calculated energy dependence of the Z-factors for p 
productiorQ, Z 1 / 1 '(E p . a), is compared in Fig.^for QGSJET- 
Ilm, E POS -LFIC, and the p arameterizations from iTan & Ng| 
(119831) and iDunerrav et alT (120031) . We consider two values, 
a = 2 and a = 3, for the slope of the CR proton spectrum, 
which bracket the physically most interesting range. 

Comparing first the results of QGSJET-IIm and EPOS- 
LHC, we observe a relatively good agreement between them 
in the high energy range for a = 2 , while somewhat larger 
differences are obtained for a steeper slope (a = 3) and for 
Ep ln < 100 GeV. This is due to the harder p production spec¬ 
tra predicted by EPOS-LHC (c.f. Figs.[[l[2] andITOl). However, 
the results for the two MC generators show a reasonable over¬ 
all agreement and, as already stated, the remaining differences 
can be used as a measure for the model uncertainties. 

Next we consider differences between the modified 
QGSJET-I I m mo del a nd the paramete r ized p spectra from 
ITan & Ngl (119831) and iDuperrav et alJ (120031) . In the en¬ 
ergy range Ep ln = 10-100GeV, where the relevant proton- 
proton interactions are covered by fixed target experiments, 
the p spectra calculated with these parameterizations agree 
approximately with those obtained using QGSJET- IIm. How- 
ever, a t hig her energies the para meterized results of lTan & Nd 
(fl983l) and lDuperrav et all (120031) are rather unreliable, which 
is best illustrated by the large difference between the Z- 


8 Here and in the following we take into account both p and n production 
when calculating Z-factors; for brevity, we use the same notation Zp P instead 
of Zp^_-. For the parameterizations bv lTan & Nd H983D and lDuperrav et af] 
120031) . we simply double the respective p yields. 
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factors calculated using the two parameter sets proposed 
by iDuperrav et al.l (120031) . Not surprisingly, the differences 
increase for larger a due to the stronger sensitivity to the for¬ 
ward p-production spectra which are less constrained by ex¬ 
perimental data. 

To get a better idea of these differences, we plot in Fig. [6] 
the ratio of /''-factors Zip' calculated with QGSJET-IIm and 
with the three parameterizations considered. One immedi¬ 
ately notices a sizeable enhancement of the p yield in the for¬ 
mer case for Ep ln < lOGeV. The differences with the results 
obtained using the parameterizations of [Tan & Ng 098 3) and 
IDuperrav et alJ (120031) reach a factor of two for Ep ln ~ 1 GeV 
and originate from the treatment of p production close to the 
kinematic threshold (c.f. Fig. ED, which is rather weakly con¬ 
strained by available experimental data. Thus, these differ¬ 
ences may be regarded as characteristic uncertainties for cal¬ 
culations of the p-yield in this low energy range. 

More importantly, for £^ ln > 100 GeV, the results 

obtain ed using the param eterizations of iTan & Nil (119831) 
and IDuperrav et al. ( 2003) start to diverge strongly from 
the results based on QGSJET-IIm, as can be clearly seen in 
Figs. [5] and [6] For instance, using the p-production spectra 
from [Tan & Nd d 19831 ), which is the standard reference for as- 
trophysical applications, one observes that the respective Zip' 
becomes practically energy-independent for El ln > 100 GeV 
(see Fig. 0, which is a consequence of the scaling picture 
used in their ansatz. However, the assumed scaling behavior 
for the inclusive particle production spectra is explicitly 
broken by the energy rise of the inelastic pp cross section 
cTpp , which has been firmly established by numerous accel¬ 
erator experiments from fixed target experiments to those 
at the LHC. Moreover, the Z-factors for p-production are 
rather sensitive to the central region in the c.m.s., because the 
p-production spectrum in pp-collisions is soft. In this region, 
the particle density rises quicker than cr ™ el , namely, as a 
power law (see, e.g., the discussion in ld’Enterria et al.1120111) . 
The resulting enhancement of p-production by high-energy 
CRs, as shown by the ratio of the respective Z-factors, is 
quite significant and may have an important impact on the 
interpretations of all sorts of CR data, not only p, by the 
PA MELA and AMS2 experiments in t his energy region (see 
e.gJBerezhko et alJ 2003HB Iasi 120091 IB Iasi & Serpico 2009 
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Cholis & Hoopeiil2014tlMertsch & Sarkarll2014l). 


Both parameterizations proposed by Duperrav et al.l (120031) 
are based on the erroneous ansatz of lLetaw et al.l (119831 ) for 
Cpp® 1 and cf'pA, i.e. they assume that the inelastic cross sec¬ 
tion is constant for E^ m > 1 GeV. The additional assumptions 
made by IDuperrav et al.l (120031) about the energy dependence 
of the very forward part of p-production spectra, which are 
responsible for the huge variations between the two param¬ 
eterizations (Fig. 0, are also questionable as we discuss in 
more detail in Appendix iBl 


5. NUCLEAR ENHANCEMENT 

Now we proceed with calculating the nuclear enhancement 
of the p-yield due to contributions of nuclear species in CRs 
and of the helium component in the ISM. As follows from 
Eq. (0 and the more detailed discussion by Kac helriefi et alJ 
(120141) . the respective partial enhancements compared to the 
yield from pp-interactions are proportional to the correspond¬ 


ing Z-factors, 


P f rr \ _ Qp i^p) _ rij Ii{Ep) Zj(E p ,ai) 

£ijlp> ~ ww) ~ ^ my z; p {E P1 a p ) ■ w 

In particular, contributions of CR nuclei may be addition¬ 
ally enhanced, if the spectral indices a* < o- 7 >, because of 
the st rong a-dependence of the Z-factors (Kachelrie fi et alJ 

l20fi . 

As follows from the general analysis dKachelriefi et al.l 
mu, in the limit of high energies and large a, one expects 
the following simple relations to hold. 


Zp(E p ,a) 

Zip {Ep, a) 
Z;\E p ,a) 
Z; p (E p ,a)~ alf 


(7) 

( 8 ) 


However, as illustrated in Fig. [T| for the case of He-p and p- 
He collisions, at finite energies and for realistic spectral slopes 
these relations are strongly modified by threshold effects re¬ 
lated to the relatively large antiproton mass. Indeed, the ra¬ 
tio Zf p /Z p p p approaches the asymptotic value 4 very slowly. 

Due to the soft p-spectrum, up to rather high energies, Zp P 
remains sensitive to p-production in the backward c.m.s. re¬ 
gion where the discussed “.4-enhancement” does not hold. 
The same effect is partly responsible for the rise of the ra¬ 
tio Z p He /Z£ p for decreasing Ep as can be seen in Fig. [7] As 
discussed in the previous Section, for small Ep ln the contri¬ 
bution of the forward c.m.s. region (moderately large zp) is 
substantially reduced by threshold effects. As a consequence, 
the relative contribution of the backward c.m.s. region Zp -C 
yj2m n/E p is increased, resulting in a “/1-enhancement” of 
the p-yield in //4-collisions. As the bottom line, we conclude 
that the nuclear enhancement of p-production is strongly en¬ 
ergy dependent in the kinematic range relevant for astrophys- 
ical applications. Hence, the use of a constant “nuclear en¬ 
hancement factor” for calculations of CR antiproton spectra 
is difficult to justify. 

In Table 0 we collect the calculated Z-factors Z l J for dif¬ 
ferent Ep ln , slopes a, and different combinations of CR and 
ISM nuclei. These results may be used for the calculation of 
secondary antiproton spectra when the partial spectra of CR 
mass groups can be approximated by a power law behavior in 
the corresponding energy range. Av an illustration , we have 
calculated the energy dependence of the nuclear enhancement 
factor for p production, e p (E p ) = JVj ef/(-Ep), with e?- de¬ 
fined in Eq. (0, in the energy range Ep ln = lOGeV-lOTeV, 
using the high energy limit of the parameterization for the 
spectra of CR species by iHonda et al. ( 2004 1: the respective 
parameters K, and a, are given in Table [2] for convenience. 
In contrast to all previous calculations, our results for e p , pre¬ 
sented in Table [3| indicate a strong energy rise of the nuclear 
enhancement for secondary antiprotons, which reaches a fac¬ 
tor of two for Ep ln ~ 1 TeV. For this particular parameteri¬ 
zation of the spectra of CR species bv \Honda et al.l ( 2004 ), the 
energy dependence of the nuclear enhancement factor ep can 
be described by a power-law, e p « 1.58 x (i?/Ge\G 0 034 . 

To understand better this result, we plot in Fig. [8] the par¬ 
tial contributions e P j(E p ) to the nuclear enhancement factor 
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from different interaction channels. One immediately notices 
a steep energy rise of the antiproton yield from interactions 
of CR helium with ISM protons. For the CR composition 
considered, the reason for this rise is three-fold. Apart from 
the trivial increase of the fraction of CR helium, the flatter 
helium spectrum compared to protons (an e < a P ) enhances 
the rise of £jj e-pi^p) because of the strong ct-dependence of 

the Z-factors Z^ (/%. a) (c.f. Eq. © and Table[2j, as noticed 
previously in iKachelriefi et al.l (120141) for 7 -ray production. 
Finally, the energy dependence of the Z-factors is affected by 
threshold effects, as we have already noticed in Fig. [7] As a 
result, the partial contributions of CR nuclei in the region of 
relatively small E p are suppressed. This suppression dimin¬ 
ishes at higher energies, where one approaches the asymptotic 
limit of Eq. ©. The latter effect dominates the energy rise 
of the contribution of heavier nuclei (CNO, Mg-Si, and Fe), 
which is shown in Fig. [S]by the dotted line. In contrast to ex¬ 
pectations, their contribution is significant in the TeV range: 
The combined antiproton yield from interactions of the heav¬ 
ier nuclei with ISM protons and helium becomes at high ener¬ 
gies comparable in magnitude with the one from the proton- 
helium channel. 


6 . CONCLUSIONS 

Accurate antiproton production cross sections are a critical 
pre-requisite for many astrophysical applications and searches 
for new physics. We have used, therefore, EPOS-LHC and 
QGSJET-II-04, two of the most advanced Monte Carlo gen¬ 
erators which reproduce numerous accelerator data includ¬ 
ing the most recent ones from LHC, to calculate the antipro¬ 
ton yield in pp-interactions. In the case of QGSJET-II-04, a 
tune of its fragmentation procedure was required for an ade¬ 
quate description of particle production at low energies and 
high thresholds. After that, we have found that the antipro¬ 
ton yields of the two QCD Monte Carlo generators agree rea¬ 
sonably well with each other and the available experimental 
data. Therefore, we conclude that the results of these two 
generators can be used to predict reliably the antiproton yield 
outside the energy range covered by fixed target accelerator 
data, E p « 10— lOOGeV. Moreover, using these Monte Carlo 
generators it is straightforward to calculate also the antipro¬ 
ton yield in pA- and .4.4-interactions. In the limiting case, 
when the spectra of CR species can be approximated by a 
power-law, we have derived the nuclear enhancement of the 
p-yield due to contributions of nuclear species in CRs and of 
the helium component in the ISM. In contrast to all previous 
calculations, our results indicate a strong rise of the nuclear 
enhancement for secondary antiprotons with energy, which 
reaches a factor of two for E p ln ~ 1 TeV. 

We have also compared our results obtained using EPOS- 
LHC and QGSJET-II-04 to the commonly used param eter- 
izations of the antiproton yield from iTan & Ng| (119831) and 


iDuperrav et al.l (120031) . In the energy range Ep ln = 10- 
lOOGeV, where the relevant proton-proton interactions are 
covered by fixed target experiments, these parameterizations 
agree approximately with the results obtained using QGSJET- 
Ilm. At higher energies, when these parameterizations are ex¬ 
trapolated outside the kinematic range constrained by experi¬ 
mental data, the differences increase fast, because the physical 
concepts used in selecting their fitting functions are incorrect. 
In particular, the assumed scaling behavior for the inclusive 
production spectra is broken by the energy rise of cr™ el , inval¬ 
idating their ansatz. The resulting increase of the p production 
should be taken into account in the interpretation of CR data, 
especially those at E > lOOGeV expected from AMS-02, 
ISS-CREAM, and CALET experiments. 

In the low energy domain E p la < lOGeV, the differences 
between the results obtained using EPOS-LHC and QG SJET- 
11-04 and the param eterizations of ITan & Ngl (119831) and 
IDuperrav et al.1 (120031) are also significant, reaching a factor 
of two at Ep ln ~ 1 GeV (c.f. Fig. [ 6 ]). The origin of these dif¬ 
ferences can be traced to the treatment of /^-production close 
to the kinematic threshold, which is rather weakly constrained 
by available experimental data. Since this energy range is the 
most interesting for dark matter searches, additional experi¬ 
mental data from, e.g., the NA61 experiment are highly de¬ 
sirable. In the absence of these data, the uncertainty of the 
antiproton yield in the low-energy range should be increased, 
reaching ~ 50% at Ef n = 1 GeV. 


The authors are grateful to Tanguy Pierog for useful dis¬ 
cussions and for his help in comparing predictions of vari¬ 
ous Monte Carlo generators to accelerator data on antiproton 
production. IVM and SSO acknowledge support from NASA 
through the grants NNX13AC47G and NNX13A092G. 


Addendum: When this paper was already prepared for 
a submission, new parameterizations for p-production in 
PiP-collisio ns have been propo s ed in two recent pub - 
lications (Idi Mauro et al.l l2014t iKappl & Winkled 120141) . 
Idi Mauro et al. (2 0141) used esse ntially the same parameter- 
ization as iDuperrav'et all (120031) and also assumed that the 
p-production spectrum scales as a power of energy (oc s A ) 
in the high energy limit. One of their proposed parameteri¬ 
zations cannot be extrapolated beyond the domain of validity 
which depends both on the energy and the slope of the pri¬ 
mary spectra. In particular, the parametrization is not suited 
for a power-law primary energy spectru m I n (E) ~ E a * 
with a p < 2.7. In turn. IKappl & Winkled (120141) followed the 
approach of ITan & Nd (11983b bv using the outdated concept 
of radial scaling, which is broken, e.g., by the energy rise of 
the inelastic cross section, as already discussed above. 


APPENDIX 

COMPARISON WITH DIFFERENTIAL P SPECTRA 

In addition to the p t -integrated spectra of p’s for pp-collisions at 158 GeV/c, plotted in Fig. [2l we present in Fig. [9] the results 
of QGSJET-IIm and EPOS-LHC for the respective differential spectra x e da/dxp / d 2 p t for fixed values of p t . Comparing the 
calculations to NA49 data (lAnticic et al.ll2010l) . we observe again a relatively good agreement, though the calculated spectra are 
somewhat harder than the measured ones. 

The same tendency is indicated by a comparison with data from the CERN Intersecting Storage Rings (ISR) (lAlbrow et al.l 
119731) at higher energies, see Fig. [TO] though any definite conclusions are hampered by both the large systematic errors of the 
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measurements and by the narrow kinematic coverage of the experiment 

Finally, the results of the two models for the momentum spectrum of p’s in pC-collisions at 12GeV/c are compared in Fig. |TT| 
to spectrometer measurements bv ISugava et al.l (119981 ). The good agreement of the spectrum calculated with QGSJET-IIm with 
the data may be regarded as an indication that the results of the model remain reasonable even when approaching the kinematic 
threshold for p-production. However, again, the narrow kinematic coverage of the experiment does not allow one to make any 
definite conclusions. 

COMPARISON OF PARAMETERIZATIONS OF P SPECTRA WITH ACCELERATOR DATA 

To trace the source of the differences between QGSJET-IIm and the parameterizations of lTan & Ng| d 19831) and lDuperrav et all 
(120031) in the predicted p-yields (Figs. [5][6]), we compare the latter with selected accelerator data in Figs. |T2l H21 and[l4] Starting 
with the results bv lTan & Nd (119831) . we observe a generally reasonable agreement with the measured p-spectra in pp-collisions 
over the wide energy range E p ~ 20 — lSOOGcVpl which dominates the Z-factors Z p ‘ for E p ' n between few GeV and few 
hundred GeV, i.e. where these factors agree approximately with those calculated using QGSJET-IIm (see Figs. [5j [6]»- However, 
the decrease of the p-yield approaching the production threshold is faster in the parameterization of lTan & Ngl ( 1983 ) than in 
QGSJET-IIm. This difference is illustrated in Fig.QAj where the respective spectra for pp-collisions at lOGeV/c are compared. 
On the other hand, in the high energy asymptotics, the scaling-like behavior of the p-spectra, implemented bvl Tan & Nd (119831) . 
is broken both by the rise of cr“ el and by the increase of the central rapidity density of secondary hadrons, as discussed in Section 
Q] It is this scaling violation which leads to a steady energy rise of the predicted Z-factors Z p p in QGSJET-IIm (c.f. Fig.0. In 
contrast, one obtains constant Z p p values using the parameterization oflTan & Nd (119831) at E p > 100 GeV. 

Let us now turn to the two parameterizations proposed by Dune rrav et al.l (120031) . Apart from neglecting the energy rise of 
the inelastic cross section, [bunerrav et al.l d2003l) made questionable assumptions concerning the energy-dependence of the spec¬ 
tral shape for p-production. In one case (parameter set 2), they assumed that the normalized (per inelastic event) p-production 
spectrum rises as a power of energy: as s A , witlfH A ~ 0 .25. While a power-law energy rise is indeed expected for hadron 
production in the central rapidity region (xp — 0) Id’Enterria et al.ll201 ll). the Feynman scaling is known to hold approximately 
for the forward spectral shape (see, e.g., a discussion in lKachelriefi & Qstapchenkoll2012l) . The outcome of this extreme assump¬ 
tion overcompensates the neglected energy rise of a™ 1 and results in a too steep increase of inclusive p-spectra, which is already 
visible in Fig. |T4] despite the narrow energy range covered by CERN ISR. In turn, in the other case (parameter set 1), the authors 
assumed a very strong scaling violation for the forward spectral shape, ex (1 — XR) constxlns (xr ~ 2E*/y/s at high energy, with 
E* being p-energy in c.m.s.). In combination with the constant cr™ el assumed, this leads to a drastic softening of the inclusive 
p-spectra in the forward direction, see Fig.fl4l 


9 The measurement has been performed with the spectrometer technique 
for a single fixed c.m.s. angle, thus covering a tiny fraction of the relevant 
kinematic space. 

1 0 Co mparing the parameterizations of p-production to the NA49 data in 
Fig- EH one has to take into account that the experimental spectrum has been 
corrected for the contributions of A and S — decays, which amounts to rescal¬ 


ing the results of lTan & Nd d 1983D and lDuperrav et afl <20031) down by some 
10 - 15% fAnticic et all20Tn) . 

11 It is noteworthy that this leads asymptotically to non-conservation of 
energy, because the energy fraction transferred to antiprotons would also rise 
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Fig. 1.— Momentum spectra of p’s in the laboratory frame in pp-collisions (left) and p Be-collisions (right) at pi a b = 19.2,GeV/c, for differen t angles 0 in 
the laborat ory frame (as indicated in the plots), calculated using QGSJET-IIm (solid, red) and EPOS-LHC (dashed, blue). The experimental data - lAllabv et alJ 
(19701) and lAmaldi et alJ (19751) . 



Fig. 2.— Feynman x-spectra of antiprotons, 1/n da/dxp (see text for definition), in pp (left) and pC (right) collisions at pi a b = 158GeV/c, calcu¬ 
lated using QGSJ ET-IIm (solid, red), EPOS-LHC (dashed, blue), and QGSJET-II-04 (dot-dashed, green), in comparison with NA49 data (Anticic et al.ll20ICt 
IBaatar et all2013l) . 
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Fig. 3.— Transverse momentum spectra of p’s at central rapidity in c.m.s. (|y| < 0.5) in pp-collisions at y/s = 900 GeV, cal culated using QGSJET-IIm (solid, 
red), EPOS-LHC (dashed, blue), and QGSJET-II-04 (dot-dashed, green), in comparison with ALICE data (Aamodt et al.l201~TI) . 
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Fig. 4.— Distribution of spectrum-weighted (a = 2.6) energy fraction Z{ (i = p, 7 ). For antiprotons, the distribution of zp = Ep/E is shown for different 
energies Ep in = 1 GeV (solid, red), 3 GeV (dashed, blue), and 10 GeV (dot-dashed, green). A similar distribution for z 7 = E^y/E for the case of 7 -ray 
production is shown for E 7 = 10 GeV by the dotted, black line marked “ 7 ”. 
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Fig. 5.— Energy dependence of Z-factors for p and n production, Zp P (Ep , a) (plotted as a function of Ep in ), for a = 2 (left) and a = 3 (right), as calculated 
with QGSJET-IIm (solid, red) and EPOS-LHC (dashed, blue), or using the parameterizations bv lTan & Nd (19831) (dot-dashed, green) and IDuperrav et all (20031) 
(dotted; the lines marked “( 1 )” (black) and “( 2 )” (red) correspond to the respective parameter sets). 
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Fig. 6 .— Ratios of Z-factors for p and n production, Z^ p , obtained using various parameterizations of p-spectra to the Z-factor calculated with QGSJET-IIm, 
for a = 2.6: solid red - ITan & Nil <19831) . dashed blue - set 1 of IDuperrav et alj <20031) , dot-dashed green - set 2 of IDuperrav et alj <20031) . 



Fig. 7.— Energy dependence of the enhancement of the Hep (solid, red) and pHe (dashed, blue) contributions to the antiproton spectrum, relative to the 
pp-case, Zp ( Ep , a)/Zp P (Ep , a) (plotted as a function of Ep in ), for a = 2.6. 
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Fig. 8 .— Energy dependence of partial contributions e^(Ep) to the nuclear enhancement factor from different interaction channels: pHe (solid, red), Hep 
(dashed, blue). He He (dot-dashed, green), and all others (dotted, black); the CR composition given in Tablets used. 




Fig. 9.— Feynman cc-spectra of p’s (c.m.s.) for fixed pt (from top to bottom: 0.1, 0.3, 0.5, 0. 7, and 0.9 GeV) in pp-collisions at pi a b = 158 GeV/c, calculated 
using QGSJET-IIm (solid, red) and EPOS-LHC (dashed, blue), in comparison with NA49 data lAnticic et al.l2010h . 
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Fig. 10.— Feynman a>spectra of p’s for fixed c.m.s. angle tan 0 = 2.66 /y/s in pp-collisions at y/s = 31 (xlO 2 ) , 45 (x 10 1 ), and 53 GeV, calculated 
using QGSJET-IIm (solid, red) and EPOS-LHC (dashed, blue), in comparison with experimental data lAlbrow et all 19731) . 
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Fig. 11.— Momentum spectra of p’s in the laboratory frame for 6 = 5.1° in p C-col lisions at pi a b = 12 GeV/c, calculated using QGSJET-IIm (solid, red) and 
EPOS-LHC (dashed, blue), in comparison with experimental data ISugava et al.l 19981) . 
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Fig. 12.— Momentum spectra of p’ s in the laboratory frame for differen t laboratory ang l es 6 (a s indicated in the plots) in pp-collisions at pi a b = 19.2 GeV/c, 
calculated using parameterizations bv lTan & Ng| 119831) (solid, r ed) or by IDuperrav et all 12003) (parameter set 1 - dashed, blue, parameter set 2 - dot-dashed, 
green), in comparison with experimental data lAllabv et alll970l : IAmaldi et al.l 19751) . 
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Fig. 13.— Feynman :c-spectra of p’s ( c.m.s.) in pp-collisions at pi a b = 158 GeV/c, calc ulate d using parameterizations bv lTan & Nd 119830 and IDuperrav et all 
120030 , in comparison with NA49 data lAnticic et all2Ol01) ; the lines are labeled as in Fig. \M 
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Fig. 1 4.— Feynman a?- spec tra of p’s for fixed c.m .s. angle tan 9 = 2.66/-v/s in pp-collisio ns at y/s = 3 1 (x 10 2 ), 45 (x 10 1 ), and 5 3 G eV: parameteriza- 
tions bv lTan & Nd d!983D and lDuperrav et all d2003D compared with experimental data dAlbrow et al!1973h : the lines are labeled as in Fig. EH 
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FIG. 15.— The pt -integrated laboratory momentum spectra of p’s in pp-collisions at = 10 GeV/c calculated using parameterizations bv lTan & Ngj 11983) 
(solid, red) and bv lDuperrav et alJ d20031) (parameter set 1 - dashed blue, parameter set 2 - dot-dashed green) in comparison with QGSJET-IIm results (dotted, 
black). 
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TABLE 1 

Z-FACTORS FOR (p + ra) PRODUCTION, (Ep, a), CALCULATED WITH QGSJET-IlM 


Ef u , GeV 

Projectile nucleus 

Target nucleus 

a = 2 

a = 2.2 

a = 2.4 

a), mb 
a = 2.6 

a = 2.8 

a = 3 

1 

P 

p 

0.0254 

0.0138 

0.00772 

0.00441 

0.00258 

0.00153 

1 

He 

p 

0.0808 

0.0442 

0.0248 

0.0143 

0.00838 

0.00501 

1 

CNO (A = 14) 

p 

0.184 

0.101 

0.0567 

0.0326 

0.0192 

0.0115 

1 

Mg-Si (A = 25) 

p 

0.273 

0.150 

0.0845 

0.0486 

0.0286 

0.0171 

1 

Fe (A = 56) 

p 

0.447 

0.245 

0.138 

0.0792 

0.0465 

0.0278 

1 

P 

He 

0.0919 

0.0498 

0.0277 

0.0158 

0.00920 

0.00546 

1 

He 

He 

0.271 

0.147 

0.0824 

0.0472 

0.0276 

0.0165 

1 

CNO (A = 14) 

He 

0.649 

0.352 

0.196 

0.112 

0.0654 

0.0389 

1 

Mg-Si (A = 25) 

He 

0.933 

0.506 

0.282 

0.161 

0.0937 

0.0556 

1 

Fe (A = 56) 

He 

1.53 

0.834 

0.468 

0.269 

0.158 

0.0944 

10 

P 

P 

0.279 

0.164 

0.100 

0.0633 

0.0413 

0.0276 

10 

He 

P 

0.979 

0.573 

0.350 

0.222 

0.144 

0.0964 

10 

CNO (A = 14) 

P 

2.67 

1.58 

0.978 

0.624 

0.410 

0.276 

10 

Mg-Si (A = 25) 

P 

4.22 

2.50 

1.54 

0.977 

0.639 

0.428 

10 

Fe (A = 56) 

P 

7.78 

4.63 

2.87 

1.84 

1.21 

0.815 

10 

P 

He 

0.970 

0.560 

0.339 

0.213 

0.138 

0.0917 

10 

He 

He 

3.16 

1.83 

1.11 

0.695 

0.449 

0.298 

10 

CNO (A = 14) 

He 

9.16 

5.33 

3.24 

2.04 

1.32 

0.875 

10 

Mg-Si (A = 25) 

He 

14.5 

8.45 

5.16 

3.26 

2.12 

1.41 

10 

Fe (A = 56) 

He 

26.0 

15.2 

9.31 

5.90 

3.85 

2.57 

100 

P 

P 

0.535 

0.308 

0.187 

0.119 

0.0789 

0.0536 

100 

He 

P 

2.03 

1.17 

0.715 

0.455 

0.300 

0.204 

100 

CNO (A = 14) 

P 

6.21 

3.61 

2.20 

1.40 

0.926 

0.628 

100 

Mg-Si (A = 25) 

P 

10.6 

6.19 

3.79 

2.42 

1.60 

1.09 

100 

Fe (A = 56) 

P 

21.8 

12.7 

7.80 

4.99 

3.30 

2.25 

100 

P 

He 

1.79 

1.02 

0.612 

0.385 

0.251 

0.169 

100 

He 

He 

6.41 

3.66 

2.21 

1.39 

0.914 

0.619 

100 

CNO (A = 14) 

He 

20.8 

11.9 

7.15 

4.50 

2.94 

1.98 

100 

Mg-Si (A = 25) 

He 

35.2 

20.2 

12.2 

7.69 

5.04 

3.40 

100 

Fe (A = 56) 

He 

72.0 

41.3 

24.9 

15.7 

10.2 

6.87 

1000 

P 

P 

0.721 

0.410 

0.248 

0.158 

0.105 

0.0715 

1000 

He 

P 

2.79 

1.60 

0.978 

0.625 

0.416 

0.286 

1000 

CNO (A = 14) 

P 

8.87 

5.12 

3.12 

1.99 

1.32 

0.903 

1000 

Mg-Si (A = 25) 

P 

15.3 

8.81 

5.37 

3.42 

2.27 

1.55 

1000 

Fe (A = 56) 

P 

31.4 

18.1 

11.1 

7.09 

4.72 

3.24 

1000 

P 

He 

2.35 

1.32 

0.787 

0.495 

0.324 

0.221 

1000 

He 

He 

8.62 

4.85 

2.90 

1.82 

1.19 

0.805 

1000 

CNO (A = 14) 

He 

29.0 

16.4 

9.81 

6.18 

4.05 

2.75 

1000 

Mg-Si (A = 25) 

He 

48.1 

27.0 

16.0 

10.0 

6.50 

4.36 

1000 

Fe (A = 56) 

He 

103 

58.1 

34.9 

22.1 

14.5 

9.85 

10000 

P 

P 

0.931 

0.516 

0.307 

0.193 

0.127 

0.0860 

10000 

He 

P 

3.60 

2.02 

1.20 

0.754 

0.493 

0.334 

10000 

CNO (A = 14) 

P 

11.5 

6.46 

3.87 

2.44 

1.61 

1.09 

10000 

Mg-Si (A = 25) 

P 

19.9 

11.1 

6.61 

4.12 

2.68 

1.80 

10000 

Fe (A = 56) 

P 

41.6 

23.4 

13.9 

8.72 

5.69 

3.84 

10000 

P 

He 

3.02 

1.65 

0.959 

0.587 

0.376 

0.249 

10000 

He 

He 

11.1 

6.14 

3.61 

2.24 

1.46 

0.983 

10000 

CNO (A = 14) 

He 

37.4 

20.5 

12.0 

7.40 

4.77 

3.20 

10000 

Mg-Si (A = 25) 

He 

63.7 

35.0 

20.5 

12.7 

8.18 

5.48 

10000 

Fe (A = 56) 

He 

135 

74.5 

43.6 

27.0 

17.4 

11.6 


TABLE 2 

Spectral parameterizations for groups of CR nuclei iIH onda et a l.I 

120041 ) 


Parameters 

H (-4=1) 

He (A=4) 

Groups of nuclei 

CNO (.4=14) Mg-Si 04=25) 

Fe (A=56) 

K 

14900 

600 

33.2 

34.2 

4.45 

a 

2.74 

2.64 

2.60 

2.79 

2.68 
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